Solitonic lattice and Yukawa forces in the rare earth orthoferrite TbFe03 
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The control of domains in ferroic devices lies at the heart of their potential for technological 
applications. Multiferroic materials offer another level of complexity as domains can be either or 
both of a ferroelectric and magnetic nature. Here we report the discovery of a novel magnetic state 
in the orthoferrite TbFeOa using neutron diffraction under an applied magnetic field. This state has 
a very long incommensurate period ranging from 340 A at 3K to 2700 A at the lowest temperatures 
and exhibits an anomalously large number of higher-order harmonics, allowing us to identify it with 
the periodic array of sharp domain walls of Tb spins separated by many lattice constants. The Tb 
domain walls interact by exchanging spin waves propagating through the Fe magnetic sublattice. 
The resulting Yukawa-like force, familiar from particle physics, has a finite range that determines 
the period of the incommensurate state. 



PACS numbers: 

Materials with magnetic transition metal and rare 
earth ions show a variety of spectacular effects originat- 
ing from the coupling between the two spin subsystems. 
The transition metal spins interact stronger and order 
at higher temperatures than the spins of rare earth ions, 
but they are also much less anisotropic. That is why 
their orientation can be controlled by the rare earth mag- 
netism. Such re-orientation transitions observed in many 
rare earth ferrites, chromites and manganites have pro- 
found effects on their magnetic, optical and elastic prop- 
erties [iHa]. 

Recently it was realized that interactions between tran- 
sition metal and rare earth spins also play an important 
role in multiferroic and magnetoelectric materials [4]-[7] . 
Thus the coupling between the Mn spins forming a spi- 
ral state in the multiferroic TbMnOa and the Ising-like 
Tb spins leads to a significant enhancement of the elec- 
tric polarization induced by the spiral |8l 19J . In GdFeOs 
orthoferrite the polarization only appears when the in- 
dependent magnetic orders of Fe and Gd sublattices are 
present simultaneously [7 , while in DyFeOs the interplay 
between the spins of Fe and Dy ions gives rise to one of 
the strongest linear magnetoelectric responses observed 
in single-phase materials [10 . 

TbFeOs is an orthorhombic perovskite (space group 
Vhnm) where Fe spins order antiferromagnetically in 
what is called G-type order along the a axis and fer- 
romagnetically (F-type) along the c axis as shown on 
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FIG. 1: Magnetic ordering of Fe^^ (brown spheres) and Tb^^ 
(blue spheres) ions in the three uniform phases of TbFeOs: 
the LT phase (panel a), the IT phase (panel b) and the HT 
phase (panel c). Also shown are the corresponding order pa- 
rameters, irreducible representations and experimental values 
of magnetic moments 11 . The various types of magnetic or- 
der depicted here are labeled as F for a ferromagnetic order- 
ing, G for a the two-sublattice antiferromagnetic Neel state, 
A for ferromagnetic ah planes stacked antiferromagnetically 
and finally C for ferromagnetic chains parallel to the c axis 
coupled antiferromagnetically. 



Fig. [T]^. This type of commensurate spin order, denoted 
as GxFz^ has an onset at approximately TN(Fe)= 650K. 
On cooling in zero field TbFeOs undergoes two transi- 
tions driven by Tb-Tb and Tb-Fe interactions JT] [T2] . 
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The ordering of Tb spins at TN(Tb)~8.5K occurs simul- 
taneously with a rotation of Fe spins in the ac plane, so 
that below 8.5K ferromagnetic components of both Fe 
and Tb spins align along the a axis, while their anti- 
ferromagnetic components are orthogonal to each other. 
The magnetic configuration of this intermediate temper- 
ature (IT) phase is F^Gz for Fe, and F!^Cy for Tb (see 
Fig.[T]3). However, below ^ 3K there is an additional spin 
re-orientation transition to a low temperature (LT) phase 
which fiips the Fe spins back to their higher temperature 
GxFz order, while the Tb spins order antiferromagneti- 
cally in the A'^Gy state (see Fig. [l^). 

Single crystals of TbFeOa were grown under 4 bar oxy- 
gen pressure using the crucible-free fioating zone method. 
Their quality was checked by X-ray diffraction. Neutron 
diffraction experiments were made on a large single crys- 
tal of TbFeOa, at the BER-II reactor of the Helmholtz 
Zentrum Berlin using the FLEX cold triple-axis spec- 
trometer with collimation of 60^-60^-60^ /c^=1.3A~^, and 
a cooled Be filter positioned in the scattered beam. Ad- 
ditional measurements were made also with the E4 two- 
axis diffractometer with A=2.8A. In both cases a mag- 
netic field was applied along the c— axis of the sample 
using a superconducting horizontal field magnet. Di- 
electric measurements were performed at the Laboratory 
for Magnetic Measurements at the Helmholtz-Zentrum 
Berlin (LaMMB), with temperatures varying between 
0.3 K and 15 K and with magnetic fields up to 1.9 T. 
Magnetic Field and temperature control was provided 
by an Oxford Instruments 14.5 T cryomagnet equipped 
with a Heliox ^He insert. An Andeen-Hagerling 2700A 
Capacitance Bridge was used to measure the capacitance 
and loss of a disc shaped sample of TbFeOa, which was 
mounted between the electrodes of a parallel-plate capac- 
itor. 

Using single crystal neutron diffraction we have probed 
the A,C,G and F-type orders in TbFeOs by tracking the 
intensity of the corresponding magnetic Bragg refiections 
in zero field and in an applied field along the c axis (see 
Methods for experimental details). In zero magnetic field 
our results are consistent with the previously observed 
sequence of the re-orientation and inverse re-orientation 
transitions [11 . Above ~8.5K we find only G-type re- 
fiections, while the development of ferromagnetic order 
is evident from the enhanced intensities of lattice Bragg 
refiections. Below 8.5K we find G- and C-type refiec- 
tions, while below 3K only A- and G-type reflections can 
be discerned. 

In an applied magnetic field (II||c) we find a far more 
complex behavior. Here we performed a series of field 
cooled measurements, while monitoring accessible A- and 
G-type refiections. In Fig. [2|a) we show in the form of a 
color plot the temperature dependence of scans along k 
around the A- type (001) refiection. At high temperatures 
this refiection is absent as there is no order of an A-type 
component for either the Fe or Tb magnetic sublattice 
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FIG. 2: Single crystal neutron diffraction data measured 
on cooling and in a magnetic field parallel to the c axis of 
H||c= 2T. All scans are measured in reciprocal space along 
(0,/c,l). In all panels except for the insert in panel (c), 
the data is plotted on a logarithmic scale in order to show 
the weaker higher harmonic reflections, (a) Temperature- 
dependent neutron diffraction measurements are represented 
in a two-dimensional plot with intensity depicted as color on 
a log scale shown on the right of the panel. White circles 
are positions of the flrst harmonic reflection computed by flt- 
ting the diffraction data at each temperature to gaussian peak 
proflles. (b) The scan measured at 3.3K is plotted here while 
the various harmonic reflections are labeled accordingly. The 
data are shown as circles while the continuous line is a flt of a 
series of Gaussian proflles to the data, (c) Scans measured at 
2.8, 2.6 and 2.1 K show here the transition from the IC phase 
as the e rapidly decreases to a smaller values into the LT-IC 
phase. At 2.8K the various higher harmonic reflections are 
clearly still evident, while at lower temperatures they merge 
closer together and appear almost as a single peak. For the 
2. IK measurement we show the flt of a single gaussian peak to 
the data which fails to account for the observed peak proflle. 
The instrumental resolution is depicted here as a horizontal 
bar. In the inset we show a flt to the same 2. IK data to a 
series of gaussians, up to the 7th harmonic, that provide a 
better model to the data. The IC wavenumber obtained from 
this flt is 6=0.002(1). 
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as indicated in Fig. 1. However, on cooling a series of 
reflections appears below 3.8K that seem to merge into a 
single peak below '^2.8K. Examination of the wavevector 
of these reflections easily establishes that they are odd 
harmonics of up to 11th order. The wavevector of the 
1st harmonic is Q = (0, e, 1) with e ~ 0.015 r.l.u. [see 
Fig.[2|b)]. The incommensurate periodicity of this phase 
that we shall refer to as IC, is approximately 67 units cells 
or -340 A. The full width at half maximum (FWHM) 
of these reflections is relatively sharp giving a coherence 
length of —700 A or approximately two full cycles of this 
unusual order. 

The physical significance of these observations is that 
the Tb-spin order in a H||c field develops a square- wave 
modulation - a periodic array of widely separated do- 
main walls in Tb magnetic order. We ascribe the ob- 
served scattering to be dominated by Tb spins first be- 
cause of its substantially greater intensity compared to 
the higher temperature Fe order found around for ex- 
ample G-type reflections and also because in zero field 
A-type reflections are associated only with Tb spin or- 
der. In the modulated A-state, Tb spins form ferromag- 
netic stripes in the ab planes with the width 170 A along 
the b axis. The a component of magnetization alternates 
from stripe to stripe, while the stacking of spins along 
the c axis is antiferromagnetic. Investigation of an Fe 
G-type reflection under the same condition suggests that 
Fe-spins are weakly perturbed by this unusual Tb-order 
(see supplementary information). 

On cooling below 2.8K Fig. |2|a) would indicate that 
the Tb modulation abruptly disappears and the Tb sub- 
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FIG. 3: The magnetic phase diagram of TbFeOs determined 
from neutron diffraction data (shown as circles) and capaci- 
tance and loss measurements (shown as triangles) that are de- 
scribed in the supplementary information (see Fig. S2-S5). In 
the insert we show the temperature dependence of the mod- 
ulation wavevector e measured in an applied magnetic field 
/jqHc = 2T. The blue line is the fit obtained using the the- 
oretical description of the IC state in terms of the periodic 
kink array. 



system returns to the zero field state with the uniform 
A^xGy order. However, closer inspection of the diffraction 
data indicates that the (001) reflection on cooling does 
not yield a simple gaussian peak shape. Rather as shown 
in Fig.[2|c), there is clear evidence of more than one peak 
at 2.8K and significant deviations from a gaussian peak 
shape at 2.1 K. This indicates that below this inverse re- 
orientation transition the state remains incommensurate 
with a much larger periodicity than above the transi- 
tion. Using carefully constrained fits of Gaussians to the 
data measured below 3. OK (see insert of Fig. [2|c)) it is 
clear that with lowering temperature e decreases sharply 
below 2.8K but does not go to zero to yield a pure com- 
mensurate (001) reflection. Rather from the temperature 
dependence of e, shown in the inset of Fig.[3j it is evident 
that below 3.7K, e decreases gradually with cooling until 
2.8K, below which it drops abruptly to a value of —0.002. 
This low temperature incommensurate (LT-IC) state has 
a periodicity of approximately 500 units cells or 2700 A, 
which distinguishes it from the higher temperature IC 
phase with a smaller periodicity. 

We have performed measurements similar to those 
shown in Fig. [2|a) at various fields as well as two zero 
field cooled measurements at 3.0 and 3.3 K where the 
field was subsequently applied isothermally in order to 
map out the various transitions that occur in this lower 
temperature regime in TbFeOa. In addition, we have 
conducted capacitance and loss measurements between 
0.3 to 10 K in a magnetic field (H||c) between and 1.9T. 
The transitions that are evident in the neutron data are 
correlated with anomalies in both the capacitance and 
loss data (see supplementary information), which allowed 
us to construct the phase diagram shown in Fig. [3| Inter- 
estingly, the capacitance and loss data imply that if the 
sample is cooled in zero field below —3 K and then field 
is applied, the transition into the LT-IC phase is not re- 
versible and this state can then be stabilized in zero field 
(see Fig. S5). 

Next we discuss the nature of interactions stabilizing 
such an unusual periodic domain wall array and holding 
the domain walls at large distances from each other. Well 
below TN(Fe)— 650K the magnitude of the ordered an- 
tiferromagnetic moment of the Fe subsystem is indepen- 
dent of temperature, while its direction in the ac plane 
described by the angle 6 can significantly vary due to the 
low magnetic anisotropy of the Fe^+ ions. In our nota- 
tions ^1 = cos^ is the order parameter of the Gx state, 
while ^2 = sin^ describes the Gz ordering. The free en- 
ergy density of the Fe subsystem is 

where the first term describes the exchange between Fe 
spins along the b axis, the second term is the magnetic 
anisotropy, which for K > favors the Gx order, and the 
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last term is the Zeeman interaction with the magnetic 
field Hz in the GxFz state. 

The free energy of Tb spins is expanded in powers of 
the order parameters 771, describing the zero phase LT 
state with antiparallel Tb spins in neighboring ah layers 
(Fig. [1^), and 772, describing the IT state with parallel 
Tb spins in neighboring layers (Fig. 
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For A = a2 — ai > the Tb subsystem would order in the 
state with 771 7^ below some temperature Tq, at which 
ai = 0. However, the interaction between the Tb and 
Fe spins favors the IT state with 7^2 7^ and ^ = ±|, 
in which both subsystems have a ferromagnetic moment 
along the a axis. Since 772 and ^2 = sin^ transform in the 
same way (see the supplemental material), this interac- 
tion is a linear coupling. 



Fe-Tb 



-A6^2. 



(3) 



For > AK, the 'unnatural' IT state with parallel Tb 
spins in neighboring layers and Fe spins rotated by 90° 
away from the easy axis, intervenes between the states 
with the 'natural' orders of Fe and Tb spins. In this way 
one obtains the zero-field phase diagram of TbFeOq pT| 

[E]. 

In addition, we consider the so-called Lifshitz invari- 
ants linear in order parameter gradients, 

/l = gi {r]idyi2 - i2dyr]i) + g2 {vidyV2 - mdym) , (4) 

which favor the experimentally observed periodic spin 
modulation along the h axis. Similar terms inducing 
modulations along the a and c axes are forbidden by sym- 
metry (see the supplemental material). Minimizing the 
total free energy - the sum of Eqs.([l][4| - we obtain the 
phase diagram shown in Fig. [4^, which includes a nar- 
row incommensurate (IC) phase region, which we assign 
to the same IC phase revealed in our neutron data. 

It is important to stress the difference between the IC 
state in TbFeOa and the long period spin spirals in non- 
centrosymmetric magnets, also described using Lifshitz 
invariants [14]. First, the crystal lattice of TbFeOs is 
centrosymmetric (inversion symmetry is only broken in 
the LT phase by the Tb spin ordering). Equation Q 
is the interaction between two distinct magnetic phases: 
the LT Tb state (odd under inversion) and the IT phase 
(even under inversion). It is only effective close to the 
boundary where these two phases have equal free ener- 
gies, which is why the IC state is observed in a very 
narrow region of the phase diagram. 

Second, spirals in non-centrosymmetric magnets result 
from the relatively weak spin-orbit coupling [14l [15] . On 



the other hand, the coupling Eq.Q likely originates from 
a stronger Heisenberg exchange: in the supplemental ma- 
terial we give symmetry arguments showing that the ex- 
change interactions between the Tb and Fe spin orders 
varying along the h axis do not cancel. Furthermore, the 
coupling between two Tb order parameters [the second 
term in Eq.Q] resulting from interactions between rare 
earth spins separated by relatively long distances, is ex- 
pected to be much weaker than interactions between the 
Tb and Fe spins described by the first term (in our cal- 
culations g2 = 0). 

Third and most important, the observation of the large 
number of Fourier harmonics in the IC state of TbFeOa 
shows that this state is qualitatively different from a mag- 
netic spiral with slowly varying spin vectors. To account 
for the difference between the isotropic Fe spins and the 
Ising-like Tb spins [12| |16j , we assumed that ci , C2 <C c 
and allowed for 25 harmonics in the periodic modulation 
of order parameters when we minimized the free energy. 
The resulting incommensurate state is shown in Fig. ^p. 
While the angle describing the Fe spins undergoes small 
amplitude fan-like oscillations around zero, correspond- 
ing to the oscillations of the weak ferromagnetic moment 
of Fe ions around the applied magnetic field H\\c^ the 
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FIG. 4: (a) Magnetic phase diagram of the Landau model 
of TbFeOs including the Fe-Tb interaction described by the 
Lifshitz invariants Eq.Q. The parameters used to obtain this 
phase diagram are: A — 0.5, K — 0.125, A — 0.275, 61 = 62 = 
1.0, bi2 2.4, gi 0.187, ^2 = 0, c = 1, and ci = C2 = 0.01. 
(b) The ^-dependence of the Tb order parameters, r]i (red 
line) and 7/2 (green line) , and the angle measured in radians 
(blue line) describing the rotation of Fe spins in the IC state 
with the period L = 340A. 
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low-temperature Tb order parameter rji exhibits sudden 
jumps. 

To understand the nature of the force that holds these 
atomically sharp domain walls at distances of ~ 170 
A from each other, we (briefly) discuss an interesting 
fleld-theoretical interpretation of the coupled system of 
rare earth and transition metal spins. Consider a sin- 
gle domain wall located at y = where the Ising-like 
LT order parameter rji shows a discontinuous jump from 
— 1 771 1 to +|?7i| or vice versa (see Fig. [5|i). Such a kink 
can be assigned the topological charge Q — (771 (+00) — 
77i(— oo))/2|77i| = ±1. The free energy per unit area of 
the domain wall is the 'bare' energy F^\^ resulting from 
interactions between Tb spins plus 



Fe = -2gO{0)Q ■ 



II' 



dy 



{K + h) 6*2 



, (5) 



where the flrst term is the Lifshitz invariant Eq.Q {g = 
gi\rii\ and ^2 = 0), describing the interaction between 
the Tb and Fe spins, while the second term is the free 
energy of Fe spins for \0\ ^ 1. Equation ([5| can be inter- 
preted as an energy of a charged plane with the surface 
charge density gQ (a 'nucleon') interacting with the fleld 
^, which describes spin waves in the Fe magnetic subsys- 
tem playing the role of pions with the mass m - 



K+h 



Minimizing Fq with respect to 6{y)^ we obtain the distor- 
tion in the Fe spin ordering produced by the Tb domain 



wall, 6{y) — 



the domain wall free energy: 



(see Fig. 
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5i), which reduces 
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When Fdw becomes negative, the domain walls 
tend to condense. Their density is, however, limited 
by Yukawa-like interactions between the domain walls, 
which result from the exchange of massive particles and 
have the range / = = '\J j^pf^- These interactions 

attract equal 'electric' charges and repel opposite ones. 
Neighboring domain walls in a periodic array have op- 
posite 'electric' charges, since topological charges of the 
domain walls alternate along the h axis. The interaction 
between two neighboring domain walls located at yi and 
y2 [see Fig. WCh)] is 



U{y2-yi) = 



\y2-y1 1 
I 



(7) 



and the total 'electrostatic' free energy of an array of 
domain walls with the charges {Qn} alternating along the 
b axis (including the 'self-energy' of the charged surfaces) 
is given by 



Fe = - ^QnU{yn -ym)Qm, 



(8) 



where yn is the position of the n-th kink. Minimizing the 
free energy density for an equidistant array of kinks (see 
Fig. |5|c)), we obtain the optimal period of the incom- 
mensurate state. Its temperature dependence flts well 
the experimental data above 2.8K, as shown in the inset 
of Fig. [3] The length scale for the period of the IC state, 
set by / ~ 150 A, is essentially the thickness of the domain 
wall in the antiferromagnetic ordering of Fe spins, even 
though such walls are not present in the IC state. Thus 
the long period of the IC state of Tb spins originates from 
the large stiffness and low magnetic anisotropy of the Fe 
magnetic subsystem. 




FIG. 5: (a) Domain wall (kink) in Tb ordering (red line) and 
the angle (blue line) describing the perturbation of Fe spins 
near the domain wall; (b) kink-antikink pair; (c) periodic ar- 
ray of domain walls with alternating charges. 

The nature of the LT-IC phase remains a puzzle, as the 
low temperature state of our theoretical model, in which 
the domain wall free energy becomes positive, is uniform. 
One possibility is that the IC order at the lowest temper- 
atures is short-ranged. Such a state might be formed by 
isolated metastable domain walls stabilized by impurities 
and crystal imperfections and separated by the average 
distance of ~1500 A. Such a conclusion can not be ex- 
cluded from the experimental data due to the small value 
of e for this state. 

The tantalizing suggestion from our work is that 
periodic domain wall arrays may be present in other 
orthoferrites and orthochromites, since they have similar 
zero fleld magnetic phases. Therefore, identifying 
domain states as the one we describe here is important 
for understanding their magnetoelectric and multiferroic 
properties. Domain walls induced by interactions 
between the rare earth and transition metal magnetic 
subsystems may also be relevant for the applications 
of this class of materials, as the dynamics of these 
domain walls, coupled by long-ranged interactions can 
strongly affect the switching of the spontaneous electric 
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polarization with an applied magnetic field and vice 
versa. 
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Supplemental material 



SYMMETRY ANALYSIS 

We use the so-called Bertaut's notations [17 to describe 
the symmetry of magnetic states of TbFeOs (see Table[l|). 
The ordering of Fe spins in the HT phase is described by 
the T4 representation. In the IT state both Fe and Tb 
spin orders have r2 symmetry, while in the LT phase the 
Fe order changes back to while the ordering of Tb 
spins has Fg symmetry (see Fig. [T]). Furthermore, the 
Lifshitz invariant has the form 



from which Eq.Q follows. We note that the rotationally 
symmetric scalar products, 



A' 



dyF 



F ■ dyA' 



and G' 



dyC 



C-dyG', (10) 



■dyF,, 



where e.g. A' • dyF = A'^ ■ dyF^ + A'y ■ dyFy + A'^ 
are also invariant under all transformations of the Pbnm 
group showing that the coupling between inhomogeneous 
rare earth and transition metal magnetic orders can orig- 
inate from Heisenberg exchange interactions [5]. 
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TABLE L Transformation properties of representations of 
Pbnm space group under the three generators of the group: 
the two glide mirrors, rhx ' {x,y,z) (1/2 — x, 1/2 + ^, 2;) 
and rhy : (x, y, z) (1/2 + x, 1/2 — 1/2 + ^), and the mirror 
ruz : {x,y,z) [x,y, 1/2 - z). 



MEASUREMENTS OF THE (Oil) G-TYPE 
REFLECTION 



Although the magnetic environment in our diffraction 
experiment constrains the portions of reciprocal space 
that we can access, we were able to also probe a limited 
region around the G-type reflection (0,1,1). In a field 
H||c= 2T we observe IC reflections below 3.5K with the 



same incommensurability e and temperature dependence 
as for the satellites around the A- type (0,0,1) reflection 
(see Fig. SI in the supplementary information). The 
commensurate (0,1,1) G-type reflection observed above 
3.5K arises only from the Fe-spin ordering (see Fig. 1). 
The fact that the intensity of the commensurate reflec- 
tion does not vary through this transition suggests that 
the Fe order is not significantly perturbed, in agreement 
with our theory. 
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FIG. SI: Single crystal neutron diffraction data measured on 
cooling and in a magnetic field parallel to the c-axis of 
H||c=2T. Scans are measured in reciprocal space along 
(0,A:,1) around the G-type reflection (0,1,1). As above the 
temperature dependent neutron diffraction measurements 
are represented in a two-dimensional plot with intensity 
depicted as color on a log scale shown on the right of the 
panel. The 1st and 3rd harmonic reflections are evident 
below 3.4 K. 



DIELECTRIC MEASUREMENTS 

Measurements of the capacitance and loss were per- 
formed both as a function of temperature between 0.3 
and 15K in field cooled mode and isothermally as a func- 
tion of magnetic field up to 1.9 T, in zero field cooled 
mode. In our measurements, a calibrated Cernox resis- 
tance thermometer (CX-1030) on the sample holder was 
used to monitor the sample temperature, which was mea- 
sured with an LakeShore 370 AC Resistance Bridge. The 
resolution of the temperature measurement is 5- 10~^, al- 
though the temperature stability decreases around 3 K, 
due to the boiling temperature of liquid ^He of 3.2 K. 
The capacitance and loss measurements were performed 
at a frequency of 1000 Hz and an excitation of 15 V. The 
resolution of the capacitance measurement is 2 • 10~^. 
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FIG. S2: Example data of capacitance and loss measured 
from a single crystal of TbFeOa between and 1.5 T as a 
function of temperature between 5 and 15 K. Here the 
sample was cooled with an applied magnetic field and 
measurements were taken continuously during warmup with 
a sweep rate of 125 mK per minute. We found that 
measurement on cooling showed exactly the same behavior. 
We also noted that that measurements after field cooling and 
zero field cooling do not differ in this temperature range. 
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FIG. S3: Example data of capacitance and loss measured 
from a single crystal of TbFeOs between and 1.9 T as a 
function of temperature between 2 and 3.5 K. Here the 
sample was cooled with an applied magnetic field while 
measurements were taken continuously during warming with 
a sweep rate of 30 mK per minute. Measurements were also 
performed down to 0.3 K, but no anomalies were found 
indicating the absence of further transitions. 
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FIG. S4: Example data of capacitance and loss measured 
from a single crystal of TbFeOa as a function of magnetic 
field between 5 and 10 K and up to 1.5 T. The sweep rate of 
the magnetic field used here was 50 mT per minute. 
Measurements were taken isothermally after zero field 
cooling. 
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FIG. S5: Example data of the capacitance and loss 
measured at 0.5 and 2.5 K up to 1.5 T. A clear peak is 
observed for data measured after zero field cooling at 0.5 T, 
consistent with the transition from the commensurate to the 
LT-IC phases observed using neutron diffraction. The same 
peak is not observed when the field was ramped down from 
1.5 T and neither is it found if we ramp the field up again. 
This behavior was consistent for a series of measurements 
between 0.3 to 2.5K and suggests that over this temperature 
range the LT-IC phase once entered into with field, is stable 
down to zero field. The sweep rate of the magnetic field for 
these measurements was 50 mT per minute. 



